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Abstract 

In this paper, we consider the gravitational radiation generated by the collision of highly relativis¬ 
tic particles with rotating Kerr black holes. We use the Sasaki-Nakamura formalism to compute the 
waveform, energy spectra and total energy radiated during this process. We show that the gravi¬ 
tational spectrum for high-energy collisions has definite characteristic universal features, which are 
independent of the spin of the colliding objects. We also discuss possible connections between these 
results and the black hole-black hole collision at the speed of light process. With these results at 
hand, we predict that during the high speed collision of a non-rotating hole with a rotating one, 
about 35% of the total energy can get converted into gravitational waves. Thus, if one is able to 
produce black holes at the Large Hadron Collider, as much as 35% of the partons’ energy should 
be emitted during the so called balding phase. This energy will be missing, since we don’t have 
gravitational wave detectors able to measure such amplitudes. The collision at the speed of light 
between one rotating black hole and a non-rotating one or two rotating black holes turns out to 
be the most efficient gravitational wave generator in the Universe. 
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I. INTRODUCTION 


In previous works we have studied the collision at the speed of light between a point 


along its symmetry axis This analyses can_c 


particle and a SchwarzscWld black hole and a point particle and a Kerr black hole 


collision between small and massive black holes [sj, the collision between stars and massive 
black holes, or the collision between highly relativistic particles like cosmic and gamma rays 


escribe several phenomena, such as, the 


colliding with black holes to name a few. We have argued in [l|, that these studies 
could give valuable quantitative answers for the collision at the speed of light between equal 
mass black holes, although we only work with perturbation theory d la Regge-Wheeler- 
Zerilli-Teukolsky-Sasaki-Nakamura, which formally does not describe this process. In fact, 
extrapolating for two equal mass, Schwarzschild black holes we obtained 1| results which 
were in very good agreement with results P| obtained through different methods. For equal 
mass black holes, we found that the collision at the speed of light between a Schwarzschild 
olack hole and a Kerr black hole along the symmetry axis gave similar results to those in 
l| . In particular, we found that 15.5% of the total energy gets converted into gravitational 
waves (for two non-rotating holes the amount is slightly less, 13%). There are as yet no 
numerical results for rotating black holes, so there is no clue as to the correctness of our 
results for rotating holes. However, should these results give radiated energies larger than 
that allowed by the cosmic censorship, and therefore by the area theorem, we would run into 
problems, and we would know that the extrapolation to equal masses is strictly forbidden. 
Our previous results do not violate the area theorem, but the mere possibility brings us to the 
study of the collision of a non-rotating black hole with a rotating one, along the equatorial 
plane, both approaching each other at nearly the speed of light. Sasaki, Nakamura and 


co-workers p have studied the infall of particles, at rest at inhnity into Kerr black holes, 
along the equatorial plane and along the axis. They found that the energy radiated for radial 
infall along the equatorial plane was about 2.6 times larger than for radial along the axis, 
for extreme holes. Accordingly, we expect the energy for highly relativistic collisions along 
the equatorial plane to be larger than that along the symmetry axis, but how larger? If the 
ratio should still be the same, 2.6, then the energy released in the collision along 

^-^axis 

the equatorial plane of a Kerr hole and a non-rotating hole at the speed of light would be 
~ 40% of the total energy. But this is slightly more than the upper bound on the efficiency 
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allowed by the area theorem, which is ~ 38.7%. Motivated mainly by this scenario, we shall 
study here the collision at nearly the speed of light between a point particle and a Kerr 
black hole, along its equatorial plane. We will then perform a boost in the Kerr black hole, 
and extrapolate for two equal mass objects, one rotating, the other non-rotating. 

Another point of interest to study this process is the hypothesis that black holes could 
be produced at the Large Hadron Collider (LHC) at CERN, which has recently been put 


Q 


forward 17| in the so called TeV-scale gravity scenarios. In such scenarios, the hierarchy 
problem is solved by postulating the existence of n extra dimensions, sub-millimeter sized, 
such that the A+n Planck scale is equal to the weak scale ~ ITeV. The Standard Model lives 
on a 4-dimensional sub-manifold, the brane, whereas gravity propagates in all dimensions. 
If TeV-scale gravity is correct, then one could manufacture black holes at the LHC. The 
ability to produce black holes would change the status of black hole collisions from a rare 
event into a human controlled one. This calls for accurate predictions of gravitational wave 
spectra and gravitational energy emitted during black hole formation from the high speed 


Q 


encounter of two particles. This investigation was hrst carried out by D’Eath and Payne 
by doing a perturbation expansion around the Aichelburg-Sexl metric, describing a boosted 
Schwarzschild black hole. Their computation was only valid for non-rotating black holes 
and it seems quite difficult to extend their methods to include for spinning black holes. Our 
approach however allows one to also study spinning black holes, which is of great importance, 
since if one forms black holes at the LHC, they will most probably be rotating ones, the 
chance for having a zero impact parameter being vanishingly small. 

Suppose therefore one can produce black holes at LHC. The hrst thing that should happen 


is a release of the hole’s hair, in a phase termed “balding” [^. The total amount of energy 
released in such a phase is not well known, and is based mostly in D’Eath and Payne’s results. 
With our results one can predict the total gravitational energy radiated in the balding phase 
when the resulting holes are rotating, a process which, as we shall see, radiates a tremendous 
amount of energy. This means that there will be a “missing” energy during the formation of 
a black hole, this energy being carried away by gravitational waves, and undetected, at least 
by any realistic present technology. Our results suggest that as much as 35% of the center 
of mass energy can be leaking away as gravitational waves, and therefore one should have 
35% missing energy. Strictly speaking, these values are only valid for a head-on collision, i.e. 
a collision with zero impact parameter, along the equatorial plane. One expects the total 
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energy to decrease if the collision is taken along a different plane, or with non-zero impact 
parameter. However, we do not expect the total energy to vary much as long as the impact 
parameter stays sufficiently small enough to give rise to black hole production. 


II. THE TEUKOLSKY AND SASAKI-NAKAMURA FORMALISM 


In this section, we give a very brief account of the Teukolsky equation, and of the Sasaki- 
Nakamura equation. Details about the Teukolsky formalism may be found in the original 
literature j^, and also in j^. For a good account of the Sasaki-Nakamura formalism we refer 
the reader to Q [ll| . 

We start from the Kerr background geometry, written in Boyer-Lindquist {t, r, 6, 0) co¬ 
ordinates: 


, 9 / 2Mr , , AMar o^,,, A sin^ d , 

ds^ = — (1- —)dt^ -— sm^ Odtdcj) + -rdr^ + Ed6^ H-— d(j) , (1) 

Zj Zj ZA Zj 

Here, M is the mass of the black hole, and a its angular momentum per unit mass. Also 
E = 4- cos^ 6 ; A = — 2Mr ] A = — Aa^ sin^ 6. 


Working with Kinnersley’s null tetrad, one can show 


that the equations for the 


Newman-Penrose quantities decouple and separate, giving rise to the Teukolsky equation. 


A-^-^{A^+^-^sR)- sVsR = -sT. 


Here, 




-K^ isKA' 2isK' 


A 


A A 

K = + a^)u! — am , 


( 2 ) 

( 3 ) 

( 4 ) 


and a prime denotes derivative with respect to r. The quantity s denotes the spin weight 
(or helicity) of the held under consideration (see for example jl^) and m is an azimuthal 
quantum number. We are interested in gravitational perturbations, which have spin-weight 
s = —2, 4-2. Solutions with s = —2 are related to those with s = 2 via the Teukolsky- 
Starobinsky identities, so we need only worry about a specihc one. For dehniteness, and 
because that was the choice adopted by Sasaki and others j^, we work with s = —2. The 
constants m, (sA depends non-linearly on a;) are separation constants arising from the 
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azimuthal function and the angular eigenfunction sZ^^{9, 0), respectively. sZ^^{9, 0) is 
the spin-weighted spheroidal harmonic, sZf^{9, 0) = sS^{9)e^'^'f’, and in turn, sSi^{9) 


satishes 
1 d 


ry^ /V 9/, 4mCOS6* , „ 

S) -I- (a^ COS^ 9 — . 2 n + ^ “I- _■ 2 /I -4 ^ ^ 0 > (^) 


. /I ^ rx\ / 9 9/1 

-T-v;^lsm0—5) + (a^ cos^^ - —^ 
sm 6* d9 d9 sm 9 


sin"^ 9 


where c = A — 2 — a^o;^ -f 2amuj. The spheroidal functions are normalized according to 
Jda\.zsi (0,0) P = 1. These functions have not been thoroughly exploited as far as we 
kn ow, but the essentials, together with tables for the eigenvalues A for s = —2 can be found 
in [3, 1A\- The source term appearing in Q may be found in j^, and depends on the 
stress-energy tensor of the perturbation under consideration. We note that, for a test particle 
falling in from inhnity with zero velocity, for example, T ~ and V is always long-range. 
This means that when we try to solve the Teukolsky eqnation 0 numerically we will run 
into problems, since usually the numerical solution is accomplished by integrating the source 
term times certain homogeneous solutions. At inhnity this integral will not be well dehned 
since the source term explodes there, and also because since the potential is long range, the 
homogeneous solution will also explode at inhnity. We can remedy this long-range nature 
of the potential and at the same time regularize the source term in the Sasaki-Nakamura 
formalism. After a set of transformations, the Teukolsky equation m may be brought to 
the Sasaki-Nakamura form [l^ : 


^X{uj,r)- T-^X{uj,r)-UX{uj,r)= C 


The functions T and lA can be fonnd in the original literature p, ll 
given by 


( 6 ) 


The source term C, is 


L = 


7oA 






and W satishes 


\Y" = _ K/Adr, 

A 


(7) 


( 8 ) 


The tortoise r* coordinate is dehned by dr^^/dr = (r^ -|- a^)/A, and ranges from — cxo at the 
horizon to -fcxo at spatial inhnity. The Sasaki-Nakamura equation (jU|) is to be solved nnder 
the “only outgoing radiation at inhnity” boundary condition (see section 4.2), meaning 


A(o^,r) = 


out^zojr* 


(X). 


(9) 


5 









The two independent polarization modes of the metric, h-|_ and hx, are given by 

.+00 


h+ +ihx = 


■'\f2TX J-c 


duj e 

l,m 






A(A + 2) - - 12a2cu2 


( 10 ) 


and the power per frequency per unit solid angle is 

J^OUt 


d?E 4a;2 ^ 

'TT < ^ 


dVtdu) TT ^ A(A + 2 ) — 12ia; — Vla?uj‘^ 


oau) 12 

‘^/ml • 


( 11 ) 


Following Nakamura and Sasaki the multipolar structure is given by 

^ (27r)V2 /■+- 

h+ + zhx =-/ U.UJ 

l,7n 

which defines implicitly. Moreover, 




( 12 ) 






(13) 


Im 


Here, _ 2 ^/m( 6 ') are spin-weighted spherical harmonics (basically 25 ^^^( 6 *)), whose properties 
are described by Newman and Penrose and Goldberg et al jl^ . 

For the befenit of comparison, the general relation between Sasaki and Nakamura’s 
X{u},r) function and Teukolsky’s R{u},r) function is given in is terms of some compli¬ 
cated expressions. One hnds however that near inhnity, which is the region we are interested 
in, this relation simplihes enormously. In fact, when r* —> cx), the relation is 


R = 


'Out^3^2a;r* 


where. 


^OUt _ 




A(A -t- 2 ) - 12ia; - I2a?uj‘^ 


(14) 


(15) 


III. THE SASAKI-NAKAMURA EQUATIONS FOR A POINT PARTICLE 
FALLING ALONG A GEODESIC OF THE KERR BLACK HOLE SPACETIME 


The dependence of the emitted gravitational wave power on the trajectory of the particle 
comes entirely from the source C in eqs. ©-CZD, which in turn through (jH)) depends on 
the source term T in Teukolsky’s equation 0 . To compute the explicit dependence, one 
would have to consider the general expression for the geodesics in the Kerr geometry Q, 
which in general gives an end result not very amenable to work with. However we hnd 
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that by working in the highly relativistic regime, the equations simplify enormously, and in 
particular for radial infall (along the symmetry axis or along the equator) it is possible to 
hnd a closed analytic and simple expression for W in (I7j)-(jHI). Accordingly, we shall in the 
following give the closed expressions for £ in these two special situations: radiall infall along 
the symmetry axis and radial infall along the equator of an highly energetic particle. 


A. The Sasaki-Nakamura equations for an highly relativistic point particle falling 
along the symmetry axis of the Kerr black hole 


We now specialize then all the previous equations to the case under study. We suppose 
the particle to be falling radially along the symmetry axis of the Kerr hole, in which case 
all the equations simplify enormously. In this case the geodesics, written in Boyer-Lindquist 
coordinates, are 


dt eo(r^ -|- a?) 
dr A 


dr 2 2 ^ 

dr ^ r'^ + 


d(f) 

dr 


A ’ 


(16) 


where the parameter cq is the energy per unit rest mass of the infalling particle. On consid¬ 
ering highly relativistic particles, en —*• oo, we hnd that the term £ takes the simple form 

in IQ) : 

/iOeoToA 


(compare with the source term 


Here, C = 
particle. 


sin^ 9 


9=0 


^ 2u;^r2(r2 + 

, the function 70 = 7 o(?") can be found in Q, 


(17) 

while /i is the mass of the 


B. The Sasaki-Nakamura equations for an highly relativistic point particle falling 
along the equatorial plane of the Kerr black hole 


We now suppose the particle is falling radially along the equator of the Kerr hole, where 
radially is dehned in the sense of Chandrasekhar 2, meaning that L = aeo, where L is 
the conserved angular momentum. In this case the geodesics, written in Boyer-Lindquist 
coordinates, are 


dt 

dr 


eo{r^ + a^) 

A 


,dr, 

'dr' 


2 2 
= ^0^ 


-A 


d(j) 

dr 




( 18 ) 
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On considering highly relativistic particles, eo —>■ C) 0 , we hnd that the source term takes 
again a very simple form (compare with the source term in ^jyi)- 

nSeo'joA 


C = 


_|_ q2^3/2 


£-*/ 


where 




f - (m - an))^ Z,„{e = h + l(™ “ 0^)1 2im(^ = h 


(19) 


( 20 ) 


2 " ' 2 ' 

Once is known, we can get the Teukolsky wavefunction near inhnity from m, and 

the waveform and energy spectra from (fTUIl -fPI). 


IV. IMPLEMENTING A NUMERICAL SOLUTION 


To compute the spin-weighted spheroidal functions, one would have to compute the eigen¬ 
values A, for each a;, using a shooting method Q or any other numerical scheme one chooses 
^ I. We have chosen to lean ourselves on work already done by Press and Teukolsky 
and to verify at each step the correctness of their results (we found no discrepancies). Press 
and Teukolsky have found, for each (/, m), a sixth-order polynomial (in auj) which approx¬ 
imates the eigenvalue A, to whithin hve decimal places (see their Table 1). This is a good 
approximation as long as au < 3, which is within our demands for I < 6. The computation 
of the spheroidal function then follows trivially, by numerically integrating 

The next step is to determine A(a;,r) from the Sasaki-Nakamura differential equation 
(jni). This is accomplished by a Green’s function technique, constructed so as to satisfy the 
usual boundary conditions, i.e., only ingoing waves at the horizon (X ~ ,r* —*• —cx)) 

and outgoing waves at in£nity(X ~ , r* —>■ cx)). We get that, near inhnity, and for infall 

along the equator (we are interested in knowing the wavefunction in this region). 


X = X'’ 


X^S 

Wr 


dr ,; X°"^ = 


X^S 

Wr 


dr,. 


( 21 ) 


The similar expression for infall along the symmetry axis is given in P . Here, X°° and X^ 
are two linearly independent homogeneous solutions of © which asymptotically behave as 


X^ ~ + B{u;)e-^^^Xr, 

—>• CX) 

(22) 

XH ^ 

► —CX) 

(23) 

Xco _ ^ 

—^ CX) 

(24) 

X“ ~ C{uj)X^^* -F D{uj)e-^‘^''Xr, 

— CX) , 

(25) 
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and Wr is the wronskian of these two solutions. Expression m can be further simplihed 
in the case of infall along the equator to 


Jj^OUt 


2iu)^B J r‘^(r‘^ + a'^yB 


(26) 


To implement a numerical solution, we hrst have to determine B{uj). We hnd B{uj) by 
solving with the right hand side set to zero, and with the starting condition = e 
imposed at a large negative value of r*. For computational purposes good accuracy is hard 
to achieve with the form so we used an asymptotic solution one order higher in l/(a;r): 


+ B(C<;) (1 - hA+2+2ama;-2ia;P) ^^-ia;r. ^ ^ _ ('27) 


Here, P = i 2 a'^uj^+i 2 m<l^\'^- 2 X^\ 2 arnuj ' These expressions all refer to the case of infall along 
the equator. For analogous expressions for infall along the symmetry axis see 

In the numerical work, we chose to adopt r as the independent variable, thereby avoiding 
the numerical inversion of r*(r). A fourth order Runge-Kutta routine started the integration 


of near the horizon, at r = r_|_ + r+e (the horizon radius r+ is r_|_ = M + (M^ 




with tipically e = 10 It then integrated out to large values of r, where one matches X^ 
extracted numerically with the asymptotic solution (|77|l . in order to hnd B{ijj). To hnd XE 
the integral in m is done by Simpson’s rule. For both routines Richardson extrapolation 
is used. 


V. NUMERICAL RESULTS AND CONCLUSIONS 

Recent studies Q, Q| on high energy collisions of point particles with black holes point 
to the existence of some characteristic features of these processes, namely: (i) the spectrum 
and waveform largely depend upon the lowest quasinormal frequency of the spacetime un¬ 
der consideration; (ii) there is a non-vanishing zero frequency limit (ZFL) for the spectra, 
and it seems to be independent of the spin of the colliding particles whereas for low- 
energy collisions the ZFL is zero; (iii) the energy radiated in each multipole has a power-law 
dependence rather than exponential for low-energy collisions. 

The present study reinforces all these aspects. In Fig. 1 we show, for an almost extreme 
Kerr hole with a = 0.999M, the energy spectra for the four lowest values of I, when the 
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FIG. 1: The energy spectra for a point particle moving at nearly the speed of light and colliding, 
along the equatorial plane, with an extreme (a=0.999M) Kerr black hole. The energy is normalized 
in units of Notice that the spectra is almost flat (for large 1 ), the ZFL is non-vanishing and 
that the quadrupole carries less than half of the total radiated energy. 


TABLE I: The zero frequency limit (ZFL) for the ten lowest radiatable multipoles. 


1 

ZFL(x^) 

1 

ZFL(x 1 ) 

2 

0.265 

7 

0.0068 

3 

0.075 

8 

0.0043 

4 

0.032 

9 

0.003 

5 

0.0166 

10 

0.0023 

6 

0.01 

11 

0.0017 


particle collides along the equatorial plane. In Fig. 2 we show the same values but for a 
collision along the symmetry axis. The existence of a non-vanishing ZFL is evident, but the 
most important in this regard is that the ZFL is exactly the same, whether the black hole is 
spinning or not, or whether the particle is falling along the equator or along the symmetry 


axis. In fact, our numerical resu 
ZFL is given by Table 1 (see also 


ts show that, up to the numerical error of about 1% the 


1| and the exact value given by Smarr 


Q). 


and this holds 


for highly relativistic particles falling along the equator (present work), along the symmetry 


Q. 


axis or falling into a Schwarzschild black hole 

The /-dependence of the energy radiated is a power-law; in fact for large I we hnd for 
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FIG. 2: The energy spectra for a point particle moving at nearly the speed of light and colliding, 
along the symmetry axis, with an extreme (a=0.999M) Kerr black hole. The energy is normalized 
in units of Notice that the spectra is almost flat (for large 1 ), the ZFL is non-vanishing and 
that the quadrupole carries less than half of the total radiated energy. 

infall along the equator 

2 2 

AE, = , a = 0.999M. (28) 

Such a power-law dependence seems to be universal for high energy collisions. Together 
with the universality of the ZFL this is one of the most important results borne out of our 
numerical studies. The exponent of I in (I2HD depends on the rotation parameter. As a 
decreases, the exponent increases monotonically, until it reaches the Schwarzschild value of 
2 {AEi ~ p) which was also found for particles falling along the symmetry axis of a Kerr 
hole. In Table 2 we show the values of the exponent, as well as the total energy radiated, 
for some values of the rotation parameter a. 


AEi = cl-^ 

a 

M 

c 

b 


0.999 

0.61 

1.666 

0.69 

0.8 

0.446 

1.856 

0.36 

0.5 

0.375 

1.88 

0.29 

0 

0.4 

2 

0.26 


Table 2. Power-law dependence of the energy radiated in each multipole Z, here shown for 
some values of a, the rotation parameter. We write AEi = p for the energy emitted for each 
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FIG. 3: The energy spectra as a function of m for I = 2 and for an highly relativistic particle 
falling along the equatorial plane of a Kerr black hole. The energy is normalized in units of 

I and AEtot for the total energy radiated away. The collision happens along the equatorial 
plane. 


This power-law dependence and our numerical results allow us to infer that the total 
energy radiated for a collision along the equator is 


AEtnt — 0.69 


M 


a = 0.999M. 


(29) 


This represents a considerable enhancement of the total radiated energy, in relation to the 


Schwarzschild case Q| or even to the infall along the symmetry axis |^. Again, the energy 
carried by the I = 2 mode {AEi =2 = 0.2^^ , a = 0.999M) is much less than the total 
radiated energy. In Fig. 3 we show, for I = 2, the energy spectra as a function of m, 
which allows us to see clearly the influence of the quasinormal modes. Indeed, one can 
see that the energy radiated for m = 2 is much higher than for m = —2, and this is due 


to the behaviour of the quasinormal frequency for different azimuthal numbers m 


Q. 


as 


emphasized by different authors [^. A peculiar aspect is that the energy radiated for m + l = 
odd is about two orders of magnitude lower than for m + I = even, and that is why the 
spectra for m = —1,1 is not plotted in Fig. 3. Notice that the ZFL is the same for m = 2 
and for m = —2 both conspiring to make the ZFL universal. In Fig. 4 we show the I = 2 
waveform for a = 0.999M as seen in 0 = 0,0 = 7r/2. By symmetry, hx = 0. 

Up to now, we dealt only with the head-on collision between a particle and a black hole. 
What can we say about collisions with non-zero impact parameter, and along the equator ? 
As is clear from our results, the fact that the quasinormal modes are excited is fundamental 
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1=2,a=0.999 



FIG. 4: The I = 2 waveform as defined by m for an highly relativistic particle falling along the 
equatorial plane of an extreme (a=0.999) Kerr black hole.The waveform is normalized in units of 
/ieo- 

in obtaining those high values for the total energy radiated. Previous studies show 

that the quasinormal modes are still strongly excited if the impact parameter is less than 

2M. We are therefore tempted to speculate that as long as the impact parameter is less 

than 2M the total energy radiated is still given by Table 2. For larger values of the impact 

parameter, one expects that total energy to decrease rapidly. On the other hand, If the 

collision is not along the equatorial plane, we do expect the total energy to decrease. For 

example, if the collision is along the symmetry axis, we know P] that the total energy is 

0.31^^ for a = 0.999M. So we expect that as the angle between the collision axis and the 

equator is varied between 0 and 7r/2 the total energy will be a monotonic function decreasing 
2 2 2 2 

from 0.69^^ to 0.31^^. Still, more work is necessary to conhrm this. Let us now consider, 
using these results, the collision at nearly the speed of light between a Schwarzschild and a 
Kerr black hole, along the equatorial plane. We have argued in previous papers m that 
the naive extrapolation fi ^ M may give sensible results, so let’s pursue that idea here. 
We obtain an efficiency of 34.5% for gravitational wave generation, a remarkable increase 
relative to the Schwarzschild-Schwarzschild collision. Now, the area theorem gives an upper 
limit of 38.7% so we may conclude with two remarks: (a) these perturbative methods pass 
the area theorem test; (b) we are facing the most energetic events in the Universe, with the 
amazing fraction of 34.5% of the rest mass being converted into gravitational waves. 
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